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Abstract 

A numerical program for the evaluation of the inclusive cross section for asso- 
ciated Higgs production with a massive weak gauge boson at hadron colliders 
is described, cr(pp/pp — > HV), V G {W, Z}. The calculation is performed in 
the framework of the Standard Model and includes next-to-next-to-leading 
order QCD as well as next-to-leading order electro-weak effects. 
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PROGRAM SUMMARY 

Program Title: vhOnnlo 
Distribution format: tar . gz 
Programming language: Fortran 77, C++. 
Computer: Personal computer. 
Operating system: Unix/Linux, Mac OS. 
RAM: A few 100 MB. 

External routines/libraries: LHAPDF (http : / /lhapdf . hepf orge . org/) , 
CUBA (http : //www . f eynarts . de/cuba/) 

Nature of problem: 

Calculation of the inclusive total cross section for associated Higgs- and W- or Z- 
boson production at hadron colliders through next-to-next-to-leading order QCD. 
Solution method: 

Numerical Monte Carlo integration. 
Running time: 

A few seconds for a single set of parameters. 



1 The program is available from 
http : //particle . uni-wuppertal . de/harlander/ sof tware/vh@nnlo 
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Figure 1: Tree-level Feynman diagrams contributing to the processes (a) qq — ¥ VH and 
(b) qq — > VHg. By crossing, one also obtains the diagrams contributing to qg —> VHq and 
qg -> VHq. 




(a) (b) (c) 

Figure 2: Sample diagrams giving rise to the 0{g 3 \ta 2 s ) terms in the cross section 



1. Introduction 

The associated production of a Higgs (H) and a weak gauge boson (V G 
{W, Z}), or "Higgs Strahlung" for short, has already proven to be an impor- 
tant process at hadron colliders (see, e.g., Refs. [U El E]). In the Standard 
Model (SM), the leading order (LO) amplitude is given by the Feynman di- 
agram shown in Fig. [I] (a). We count it as order g 4 , where g is the weak 
gauge coupling. Through next-to-leading order (NLO) QCD, 0(g 4 a s ), it can 
be represented as the convolution of the production of a virtual gauge boson 
V* with the decay rate of V* into the on-shell gauge boson V and the Higgs 
boson H, schematically: oyu = &V* <8> ^v*-^vh- The precise formula will be 
given below. 

The bulk of the next-to-next-to-leading order (NNLO) QCD corrections, 
0{g 4L a 2 s ) ) is given by the same formula and will be referred to as DY (for 
"Drell-Yan like") terms in the following. 
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Figure 4: Typical diagram for the ggHZ contribution. 
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At 0(a 2 8 ), there are contributions to the cross section which cannot be 
derived from the Drell-Yan process as before. One such class is given by 
diagrams where the Higgs boson is radiated off a virtual top-quark loop 
which is inserted either in a real or a virtual gluon line. Sample diagrams are 
shown in Fig. [2] Apart from that, for ZH production, both the Z and the 
H can couple to the top-quark loop; examples are shown in Fig. [3] All these 
diagrams enter the total cross section at 0(g 3 \ t al) through the interference 
with the LO amplitudes shown in Fig.fTJ Although these are not the only 
contributions induced by top-quark loops, we will refer specifically to these 
0(g 3 \ t a 2 ) terms as the "top" terms. 

For ZH production, yet another top-quark induced contribution enters at 
0(a 2 ). Sample diagrams are shown in Fig.[4j They differ from the previous 
ones in that they contain two initial state gluons. Since this does not occur at 
lower orders in a s , it is the square of the corresponding amplitude that enters 
the cross section, rather than the interference with lower order terms. It will 
be referred to as the "ggHZ" contribution in what follows. Note that there 
are also DY diagrams with gg initial state, but they belong to the double real 
emission terms (e.g. gg — > VHqq), while the ggHZ contribution is a purely 
virtual correction. In our counting of coupling constants, it enters the cross 
section at 0(g 2 X 2 a 2 ). 

The NLO QCD corrections have been known for a long time [3]. The 
DY NNLO corrections, based on the NNLO results for the actual Drell-Yan 
process [HIE] and including the ggHZ terms [7], were first presented in Ref. [8]. 
The top effects were subsequently evaluated in Ref. [9]. For a Higgs mass of 
Mh = 125 GeV, inclusion of the NLO corrections enhances the total inclusive 
cross section for WH production by about 31% (27%, 41%) at the LHC8 
(LHC14, Tevatron); the effect of the NNLO DY terms is an additional 3% 
(3%, 10%), and the top effects amount to about 1% in all three cases. For ZH 
production, these numbers are very similar; the additional ggHZ component 
amounts to about 5% (9%, 0.2%). 

These numbers show that the numerical effect of the NNLO DY corrections 
is rather small; however, including them significantly reduces the theoretical 
uncertainty due to scale variations [8] , in particular for the WH process. 
For ZH production, on the other hand, the fact that the ggHZ terms are 
separately finite, gauge independent, and proportional to a 2 (fi) introduces 
an additional and significant amount of theory uncertainty (see also Refs. [T(H 

HI]). 

Currently, the program vhOnnlo includes the DY, the top, and the ggHZ 
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contributions at 0{a 2 s ). Electro-weak corrections to the Higgs Strahlung 
process are of the order of 5% and thus larger than the QCD scale uncer- 
tainty [T2l H3] . Since they do not depend on any of the QCD parameters like 
PDFs or the strong coupling, we included them in vh@nnlo in the form of a 
correction factor by linearly interpolating the numbers given in Ref. [H] . 

Note that vh@nnlo only evaluates the total inclusive cross section; kine- 
matical distributions, for example in the Higgs' transverse momentum or 
rapidity cannot be calculated with this program. We refer the reader to 
Refs. JT5l [IH [T7] (and references therein) for that purpose. 

2. Description of the program 

2. 1 . Structure 

This section describes the internal structure of vh@nnlo. 
2.1.1. DY terms 

The core of the program consists of the Fortran 77 code zwprod by 
W. van Neerven which is publicly available from Ref. [IS]. It evaluates the 
production of an off-shell weak gauge boson V* of invariant mass q 2 through 
NNLO QCD by implementing the results for the partonic cross section <r^g 
for the process a(3 — » V* + X (a, (3 G {q, q, g}) of Refs. [5j [6] and integrating 
them over the parton density functions (p a / p : 



The DY terms of the inclusive total cross section for VH production are then 
obtained by folding this expression with the decay rate of the off-shell vector 
boson dr: 




(1) 




(M H +M V ) 2 



dg 2 ov(g 2 ,s) 



dT(MlM 2 ;q 2 ) 
dq 2 



(2) 



where 



dT(M 2 ,M 2 ;q 2 ) _ G F M± X^jM 2 , M 2 ; q 2 ) 
dq 2 " 2V2tt 2 (q 2 - M 2 f 



1 + 



A(My, M^; q 2 ) 
M 2 /q 2 




(3) 
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and 



A( W )^(l---^) 2 -4^. (4) 

V z z/ z 

In vh@nnlo, the program zwprod is modified in such a way that it per- 
forms the three integrations of Eq. ([TJ and ^ simultaneously using the Vegas 
algorithm [19] as implemented in the CUBA library |20j . 

For the PDFs, vh@nnlo links the LHAPDF library [2T] which allows one to 
switch between different PDF sets conveniently. Neither LHAPDF nor CUBA 
are part of the distribution of vhOnnlo; they have to be installed separately. 

2.1.2. ggHZ terms 

The ggHZ terms of order 0(g 2 X 2 a 2 ) are calculated by a routine that was 
generated with the help of FeynArts [22] . Due to the top-quark loop, one 
arrives at massive box integrals that are evaluated using a slightly modi- 
fied version of the LoopTools/FF library (23J El] which is included in the 
distribution of vh@nnlo. For the PDFs, again LHAPDF [21J is used. 

2.1.3. top terms 

The top contributions of order 0(g 3 X t a 2 ) are evaluated by a C++ routine. 
The calculation is based on the heavy-top limit which was shown to work 
well within the overall available accuracy [9] . Folding the partonic differential 
cross section with the PDFs and integrating it over the phase space requires 
a three (six) dimensional numerical integration for the virtual (real) contri- 
bution, which is again performed using LHAPDF and the CUBA implementation 
of VEGAS. 

2.1.4- Electro-weak corrections 

As already mentioned above, electro-weak corrections are implemented in 
the form of a two-dimensional grid spanning over the values Mr- G [80, 200] GeV 
in steps of 5 GeV, M H G [200, 300] GeV in steps of 10 GeV, and G 
{7, 8, 9, 10, 14} TeV. The correction factor is very flat; at intermediate val- 
ues it is therefore obtained by linear interpolation. 
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2.1.5. Total cross section 

In summary, the individual orders of the inclusive Higgs Strahlung cross 
section are calculated as follows: 

o"lo = (1 + 5ew)o"lo , 
o"nlo = (1 + ^ewVnlo > (5) 

^NNLO = (1 + (WVnNLO + ^ ggHZ + ^ tOP , 

where <5ew denotes the electro-weak correction factor. 

2.2. Installation and compilation 
2.2.1. Installation 

Unpacking the tar ball will create the following directory structure: 

vMnnlo/ 
inf iles/ 
mainf iles/ 
manual/ 
output/ 
source/ 
x/ 

vh@nnlo: this will also be referred to as the root directory of vh@nnlo. It 
contains the file README which briefly describes the most important top- 
ics for installation, compilation, and operation of the program vh@nnlo. 

source : contains both the sources of vh@nnlo, as well as the main Makefile 
and a very simple compilation script named vh@nnlo_install. 

x: should be empty upon first installation. It will contain the executable 
once vhOnnlo is compiled. 

output : a directory to store output files. 

inf iles: a directory to store input files. 

manual : contains this manual. 

mainf iles: a directory to store the main program. Currently, it contains 
only the file main.f. 
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2.2.2. Compilation 

Before compiling, it might be necessary to adjust some of the library 
paths in source/Makefile. 

Full compilation is most easily done by saying 

. /vh@nnlo_install 

in the directory source. This will compile the LoopTools/FF library, and 
the ggHZ, DY, and top components of the calculation. The main executable 
x . main will be copied to the directory x and should be called from the root 
directory of vh@nnlo. 

Sometimes it is obvious that the ggHZ component does not contribute to 
the cross section. This is true for WH production, for example, or if only 
the NLO result is needed. In that case, one may work with only a subset 
of the code and exclude, for example, the LoopTools/FF library from the 
compilation. Instead of running the installation script, one may then simply 
say 

make GGHZ=no 

in the source directory. 

2. 3. Operation 
2.3.1. Input 

vhOnnlo is controlled by a single input file. An example named in. init 
is included in the distribution in the directory inf iles. It is structured like 
an SLHA input file [25J, but is not SLHA compatible (it uses its own Blocks, 
for example). The content of in. init as contained in the distribution is: 
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Block REAL 

1 8.000000D+03 # sqrt(s) [GeV] 

2 125.0000D+00 # mh [GeV] 

3 0. 100000D+01 # muR/q 

4 0.100000D+01 # muF/q 

8 172. 5d0 # mt [GeV] 

9 4.75d0 # mb [GeV] 

11 0.911876D+02 # Mz [GeV] 

12 0.803980D+02 # Mw [GeV] 

13 0.249520D+01 # GammaZ [GeV] 

14 0.214100D+01 # GammaW [GeV] 

15 0.116637D-04 # GFermi [l/GeV~2] 
18 0.508000D-01 # sin~2(thetaC) 

Block INTEGER 

1 1 # [l:pp] [2:ppbar] 

2 2 # [0:L0] [1:NL0] [2:NNL0] (adjust PDF set as well!) 

3 1 # [0:WH] [1:ZH] 

4 # PDF number 

5 5 # nf 

6 1 # [0:w/o] [l:with] electro-weak corrections 
Block CHAR 

1 MSTW2008nnlo68cl . LHgrid # PDF set 

It is divided into "Blocks" that are characterized by the type of the input 
variables they contain (REAL*8, INTEGER, and CHARACTER). Each line inside 
a Block consists of (from left to right) 

• at least one whitespace character 

• a label of type INTEGER 

• the value of the input variable 

• (optionally) a comment, introduced by a hash symbol (#) 

Quite generally, vh@nnlo will ignore any text to the right of the hash sym- 
bol #. 

In the following, the individual input parameters are described in more 
detail. All masses, energies, and decay widths are to be given in GeV. 
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Block REAL — input variables of type REAL*8 

1 y/s, the hadronic center- of-mass energy 

2 — Mh, the mass of the Higgs boson 

3 Pr/ \f^t ''> the renormalization scale relative to a/ (pv + Ph) 2 
Pv and are the 4-momenta of the final state weak gauge 
boson and the Higgs, respectively. At LO, q 2 = (p v + Ph) 2 is 
thus the square of the partonic center-of-mass energy. 

4 ^f/^/q 2 , the factorization scale relative to y/(pv + Ph) 2 
See entry 3. 

8 — M t , the on-shell top quark mass 

9 — Mb, the on-shell bottom quark mass 

11 — Mz, the Z boson mass 

12 — Mw, the W boson mass 

13 — r z , the Z boson decay width 

14 — Y w , the W boson decay width 

15 — Gp, the Fermi constant (in GeV~ 2 ) 

Note that the electro-weak couplings g, g' are expressed in 
terms of Gp, Mw, and Mz; in particular, we use Gf = 
7ra QED /(v / 2M^ sin 2 6 W ) and sin 2 9 W = 1 - M w /M|. 
18 — sin 2 9 C , the Cabbibo angle 

Mixing of the first and second quark generation with the 
third generation is neglected, i.e. vh@nnlo assumes V qt = S q b, 
q G {d,s,b}. The error introduced by this approximation is 
completely negligible. 
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Block INTEGER — input variables of type INTEGER 

1 — collider type: 1 = pp, 2 = pp 

2 order of calculation: = LO, 1 = NLO, 2 = NNLO 

This parameter affects only the order of a s that is taken into 
account for the partonic cross section. It has no influence 
on the order at which a s or the PDFs are evolved. This 
is determined by the choice of the PDF set (Block CHAR, 
entry 1). 

3 — type of process: = WH, 1 = ZH 

4 — number of the PDF set member 

Modern PDF sets usually provide a means to evaluate the 
corresponding uncertainties. For this, they contain several 
"member" PDFs within one set. The set itself is specified in 
Block CHAR, entry 1, see below. 

5 number of flavors taking into account in the running of a s 

6 include electro-weak correction factor: = no, 1 = yes 

Block CHAR — input variables of type CHARACTER 

1 — the name of the PDF set as defined in LHAPDF 

The perturbative order of the PDF set will determine the 
order of the DGLAP evolution of the PDFs, and also the order 
to which a s is evolved from Mi to /ir. 

In order to run the program with this input file, say 

x/x.main inf iles/in. init 

in the main directory of vh@milo. 

2.3.2. Output 

The input file of section 2.3.1 will lead to the following output, written to 
the file output /out .vh: 
# 

# +++ vhOnnlo version 1.20 +++ 

# Authors: T. Zirke and R.V. Har lander (Drell-Yan part) 

# 0. Brein (ggHZ part) 

# Based on zwprod.f by W. van Neerven and the paper 

# 0. Brein, A. Djouadi, R.V. Harlander, 

# Phys.Lett. B579 (2004) 149, hep-ph/0307206 . 

# Uses LoopTools by T. Hahn and FF by J. van Oldenborgh. 
# 
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# In addition, please cite the following papers: 

# Han:1991ia 

# Hamberg: 1991np 

# Harlander : 2002wh 

# Brein:2003wg 

# Ciccolini:2003jy 

# Kniehl:1990iv 

# 

Block REAL 



1 


0. 


800000D+04 


# 


sqrt(s) [GeV] 


2 


0. 


. 125000D+03 


# 


mh [GeV] 


3 


0. 


. 100000D+01 


# 


muR/q 


4 


0. 


. 100000D+01 


# 


muF/q 


8 


0. 


, 172500D+03 


# 


Mt [GeV] 


9 


0. 


475000D+01 


# 


Mb [GeV] 


11 


0. 


.911876D+02 


# 


Mz [GeV] 


12 


0. 


803980D+02 


# 


Mw [GeV] 


13 


0. 


249520D+01 


# 


GammaZ [GeV] 


14 


0, 


214100D+01 


# 


GammaW [GeV] 


15 


0, 


. 116637D-04 


# 


GFermi [GeV] 


16 


0. 


.755625D-02 


# 


alphaQED 


17 


0. 


222646D+00 


# 


sin~2(thetaW) 


18 


0. 


508000D-01 


# 


sin~2(thetaC) 



Block INTEGER 

1 1 # [l:pp] [2:ppbar] 

2 2# [0:L0] [1:NL0] [2:NNL0] 

3 1 # [0:WH] [1:ZH] 

4 # PDF number 

5 5 # nf 

6 1 # [0:w/o] [1 :with] electro-weak corrections 
Block CHAR 

1 MSTW2008nnlo68cl.LHgrid # PDF set 
Block SIGMA 

1 0.398313D+00 # sigma(all) [pb] 

10 0.117070D+00 # alpha_s(mz) 

11 0.399203D+00 # sigma(DY) [pb] 

12 0.154259D-01 # sigma(gg->HZ) [pb] 

13 -.510000D+01 # deltaEW ['/,] 

14 0.404315D-02 # sigma(top) [pb] 
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One observes that vh@nnlo includes the input to the output file, and adds 
a new Block, named SIGMA. Its entries are given as follows: 



Block SIGMA — output o/vh@nnlo 

1 o~tot, the total inclusive cross section, including all corrections 

requested in the input file 

10 a s (Mz), the input value for the strong coupling constant, de- 
fined by the PDF set 

11 — cr DY , the DY part of the cross section 

12 — o~ ggHZ , the ggHZ component of the cross section 

13 — Sew, the electro-weak correction factor 

14 — a top , the top component of the cross section 

Note that vh@nnlo will output the results only up to the order requested 
in the input file. Thus, for example, cx ggHZ and o~ top (entries 12 and 14) will 
be zero unless the requested order of the calculation is NNLO (i.e. Block 
INTEGER, entry 2 is equal to 2). 

2.3.3. Citations 

vhOnnlo collects a number of results which need to be acknowledged in a 
proper way. It is insufficient for the user of vh@nnlo to simply refer to this 
manual without explicitly referencing the relevant literature that entered 
into the calculation. In order to facilitate this task, vh@nnlo outputs a list 
of "texkeys" in the header of the output file which have to be referenced (see 



Section 2.3.2). The full references can be obtained with these texkeys from 



the INSPIRE database using the find texkey command. 
2.3.4- Parameter scans 



The input file described in section 2.3.1 defines one single set of parame 



ters. In order to perform parameter scans, one could modify the main pro- 
gram mainf iles/main.f . However, we recommend leaving the Fortran code 
of vhOnnlo untouched, and rather operate with scripts such as slharoutines 
which are available from Ref. 



3. Conclusions 

vhOnnlo is a program that collects the most up-to-date results for inclu- 
sive production of Higgs bosons in association with a weak gauge boson in 
the SM. Individual perturbative contributions can be separately considered. 
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Various parameters can be adjusted. In particular, parton densities can be 
easily changed through the use of the LHAPDF library. 

We hope that the program will be useful for Higgs physics at the LHC. 

Acknowledgments. This work was supported by DFG, project HA 2990/5-1, 
and by the Helmholtz Alliance "Physics at the Terascale" . 
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